Non-Abelian SU (2) gauge fields through density-wave order and strain in graphene 
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Spatially varying strain patterns can qualitatively alter the electronic properties of graphene, 
acting as effective valley-dependent magnetic fields and giving rise to pseudo-Landau-level (PLL) 
quantization [1]. Here, we show that the strain-induced magnetic field is one component of an 
SU(2) non-Abelian gauge field within the low-energy theory of graphene, and identify the other two 
components as period-3 charge-density waves. We show that these density- waves, if spatially varied, 
give rise to PLL quantization. We also argue that strain-induced magnetic fields can induce density- 
wave order in graphene, thus dynamically gapping out the lowest PLL; moreover, the ordering should 
generically be accompanied by dislocations. We discuss experimental signatures of these effects. 



The discovery of graphene [1] — a carbon monolayer 
with low-energy electronic properties governed by the 
Dirac equation [2] — has stimulated enormous interest 
in condensed matter systems having Dirac quasiparti- 
cles. Although other systems supporting Dirac quasi- 
particles have subsequently been discovered (e.g., the 
surface bands of topological insulators [3, 4]), graphene 
is uniquely tunable through lattice deformations and 
strains, being a soft two-dimensional membrane [2]. 
Strains alter the electronic structure of graphene by mod- 
ulating the hopping amplitudes between neighboring lat- 
tice sites. A striking example of the electronic conse- 
quences of strain is that certain spatially varying pat- 
terns of deformation can mimic the effects of a pseudo- 
"magnetic field" that has opposite signs in the vicinity 
of the two Dirac points (which we shall refer to as the 
valleys K and K' [5]). Such pseudo-magnetic fields, on 
the order of 300 Tesla, were recently demonstrated in pi- 
oneering experiments [6] on nanoscale graphene bubbles; 
more recently, effective fields of order 60 Tesla were real- 
ized by artificially designing a lattice of carbon-monoxide 
molecules on copper [7]. At these fields, each valley is 
deep in the quantum Hall regime, so that its electronic 
structure consists of well-spaced pseudo-Landau levels 
(PLLs). Because the PLLs are highly degenerate, one 
expects correlation effects to be strong within them; in- 
deed, recent works have shown that, in the presence of 
interactions, partially-filled PLLs are unstable to form- 
ing ordered states such as valley ferromagnets, spin-Hall 
phases and triplet superconductors [8, 9]. 

In the present work we show that the strain-induced, 
valley-dependent magnetic field is one component of a 
non-Abelian SU(2) gauge field within the low-energy the- 
ory of graphene. We identify the other two generators of 
this SU(2) gauge field as period-3 charge-density waves 
(3CDWs) (Fig. 1) that mix the K and K' valleys. We 
show that these charge-density waves act as gauge poten- 
tials: when their amplitude is constant, they move the 
Dirac cones [Fig. 2(B)]; but when their amplitude is spa- 
tially varied, they can give rise to Landau-level quantiza- 
tion, as shown in Fig. 3. Although methods for realizing 
non-Abelian SU(2) gauge fields had previously been pro- 



posed in ultracold atomic settings [10-12] and in twisted 
bilayer graphene [13], such fields have yet to be experi- 
mentally realized in condensed matter. The present work 
suggests an alternative approach, which might be easier 
to implement, e.g., in molecular graphene [7]. 

Having established the SU{2) gauge structure, we turn 
to the effects of the 3CDW patterns on the PLL struc- 
ture induced by strain. We find that, although these per- 
turbations do not open up gaps in unstrained graphene, 
they do gap out the lowest (i.e., zero-energy) PLL. On 
general grounds, then, we expect these gaps, and the 
corresponding 3CDW patterns, to be dynamically gener- 
ated by electron-electron or electron-phonon interactions 
whenever graphene is strained, as they would reduce 
the ground-state energy. (The relation between different 
mass gaps and the corresponding ordered states was pre- 
viously explored, for unstrained graphene, in Refs. [14- 
16].) 

After discussing the 3CDW patterns, we turn to their 
defects (i.e., dislocations), and show that these defects 
are entwined with the ordering in a distinctive way, owing 
to the valley-dependence of the pseudo- magnetic field. 
In contrast with the case of a regular field, for a strain- 
induced field a uniform 3CDW perturbation does not mix 
the spatially coincident Landau orbitals in the two val- 
leys, as these are counter-propagating. However, a 3CDW 
perturbation with a defect at the origin can mix the val- 
leys and open up a gap. Thus, in experimental geometries 
such as that of Rcf. [6] , the defects as well as the order 
are likely to be dynamically generated. 

Model. In the absence of interactions, the tight-binding 
Hamiltonian of strained graphene reads 

Ho = J2 J2 {t + St a (T i ))(a\r i )b(T i + S a ) + h.c), (1) 

r 4 a=l,2,3 

where <^ ( r i) is the strain-induced variation of the 
nearest neighbour hopping amplitude between the A- 
sublattice site at and the i?-sublattice site at + d a of 
the bipartite honeycomb lattice [2]. The vectors S a con- 
nect any A-sublattice atom to its three B-sublattice near- 
est neighbors. In the absence of strain, the low-energy 
excitations correspond to linearly dispersing states close 
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FIG. 1. 5(7(2) gauge potentials [i.e., charge-density-wave 
(3CDW) patterns] and their defects. Panels (a) and (b) show 
the patterns corresponding to £ Xi i and C, y> i respectively (see 
Table I for definitions); the other two patterns are related 
to these by translation. Black (white) circles denote excess 
positive (negative) charge; larger circles correspond to greater 
excess charge. Panel (c) sketches a dislocation in the 3CDW 
pattern; the dislocation core is marked in gray. Dashed red 
lines trace the crests of the 3CDW pattern. 



to the two Dirac points at momenta ±K with K = 
(47r/3-\/3ao) e x, «o being the carbon-carbon bond length 
[2]. Near the Dirac points K and K' the wavefunctions 
of such states can be written as four-component spinors 
"F = (Vm.,K) V'b.Kj Vu.K'i V'-BjK') where the first index de- 
notes the component of the wavefunction on the A(B) 
sublattice of the honeycomb unit cell, and the second 
index denotes the component of the state that is asso- 
ciated with the K (K') valley. The low energy effective 
Hamiltonian close to the Dirac points reads as: 



H = v F [p x T x + PyTy] 



(2) 



where T x — t 3 (Ti, T y = tq<72 , vp is the Fermi velocity, 
and the a and r operators are Pauli matrices acting on 
sublattice and valley indices respectively. We have not 
included the physical spin index as it does not affect our 
analysis, provided that the spin-orbit coupling is negligi- 
ble. 

There are three terms in the low-energy theory (i.e., 
"charges") that commute with the Hamiltonian (2): 
Qi = -T2<72, Q2 = T~ia 2 , Q3 = T 3 a - These realize 
an SU(2) pseudo-spin algebra [Qi,Qj] — 2ieijkQj- We 
also define the electromagnetic U(l) charge Qo = t <t 
(i.e., the identity operator), which commutes with the 
other charges. We can minimally couple Ho to the gauge 
potentials associated with these charges, thus arriving at 
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TABLE I. Density-wave-based gauge potentials in graphene. 
This table lists the microscopic and low-energy forms of the 
perturbations, as well as the forms after projection onto the 
lowest pseudo-Landau level (LPLL). 



the Hamiltonian: 
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We turn to the microscopic origins of the SU(2) gauge 
potentials, Cti = T^Q iy where fj, — x,y and i = 1,2,3. 
Of these, £^,3 comprise the familiar strain-induced vec- 
tor potential. Strain generates a gauge field given by 
A x + iA y = Y, a =i 2 3 St a (r)e ±iK '° a near the Dirac points 

±K. Note that J2 a =i 2 3 ^a(r)e ±jK ^ a is complex be- 
cause the nearest-neighbor hoppings are not symmetric 
under inversion. The real part of the strain gauge field 
A° is the same in both valleys and therefore couples to 
Qo! it can be gauged away assuming time-reversal sym- 
metry holds. On the other hand, the imaginary part iA y 
has opposite sign in the two valleys and couples to Q3 
leading to the valley-dependent magnetic fields realized 
in the experiments of Ref. [6] . 

The four remaining gauge potentials (see Table I) orig- 
inate, as we shall now see, as 3CDWs. That they should 
be charge density waves can be seen as follows: (a) the 
perturbations mix the valleys, and must therefore involve 
spatial modulations that enlarge the unit cell; (b) they 
do not mix the sublattices (i.e., they are proportional 
either to 03 or to <tq), and can therefore include only on- 
site charge offsets and intra-sublattice (e.g., next-nearest 
neighbor) hopping. Two simple perturbations satisfying 
both criteria are charge modulations of wavevector G 
(where G = K — K' is a vector connecting the two Dirac 
points), which realize t\ and ti respectively: 



Ti <H> cos(G • r), T2 f> sin(G • r) 



(4) 



The gauge potentials (,x,i are realized when the density 
waves on A and B sublattices are it out of phase, whereas 
the potentials are realized when the density waves on 
the A and B sublattices are in phase. Fig. 1 shows the 
two corresponding density- wave arrangements, which are 
listed in Table I. 

Numerical band structure calculations on nano- ribbons 
including these terms are shown in Figs. 2 and 3; as one 
can see, the £'s do not open up gaps in the absence of 
strain, but do shift the Dirac points in momentum space, 
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FIG. 2. (A) The band structure in the presence of strain, cal- 
culated for a 100-site thick graphene nanoribbon with zigzag 
edges, showing pseudo-Landau levels (PLLs). All energies 
are computed in terms of the unperturbed hopping matrix 
element; in the presence of strain, the hopping along one of 
the three directions increases linearly from t to 1.9t across 
the sample. Note that the zero-energy level is in fact four- 
fold degenerate, including both two zero-energy PLLs and 
two modes localized on zigzag edges. (B) The location of 
the Dirac cones in the absence of strain, both with the per- 
turbation £ Xii (thin solid line) and without it (thick line). 
As noted in the main text, constant gauge potentials are ex- 
pected to shift the Dirac cones. (C) The band structure in the 
presence of both strain and the perturbation C, x ,i- Note that 
the fourfold-degenerate zero-energy mode is completely split, 
with the two more widely separated levels corresponding to 
the lowest PLL. (This was checked by explicitly computing 
the wavef unctions.) 



as a vector potential is expected to do. Moreover, if the 
coefficient of Cy^i say, is varied linearly with x (cf. the 
Landau gauge description of a uniform magnetic field), 
it gives rise to PLL quantization as shown in Fig. 3. We 
have thus established that these terms can be regarded, 
together with the strain-induced gauge potentials, as en- 
abling the realization of a general SU(2) gauge field. 

Strain-induced PLLs and their mass gaps. The pseu- 
dospin 5/7(2) symmetry (at low energies) allowed us to 
treat the strain and 3CDWs on the same footing in the 
preceding discussion. We now break this symmetry by 
considering the PLL structure created by a strain pattern 
inducing a uniform pseudo-magnetic field. The eigen- 



FIG. 3. Pseudo-Landau levels (PLLs) generated by the SU(2) 
pseudo-gauge potential A x Qi = BQiy on a 100-site nanorib- 
bon with zigzag edges. The value of Q2 varies linearly from 
zero at one end of the ribbon to 0.74i at the other end. The 
main panel shows wavefunctions in the first, second, and third 
distinct PLL, calculated on a 100-site nanoribbon with zigzag 
edges. The inset shows the overall level structure (cf. Fig. 2A, 
for the case of regular strain), demonstrating that the PLLs 
are indeed flat for momenta near the Dirac cone. 



states of T-L then fall into PLLs at energies E n <~ y/n, 
and in particular there is a PLL at zero energy in each 
valley [2], which we term the LPLL. While the LPLL 
shares some features with the zero-energy Landau lev- 
els induced by a real magnetic field, it is distinct in 
two essential ways, as follows. (1) In contrast with the 
case of a real magnetic field (in which the two lowest 
Landau level wavefunctions in the two valleys are lo- 
cated on opposite sublatticcs), the LPLL wavefunctions 
in both valleys are located entirely on the A sublattice [8] . 
(2) The wavefunctions in the LPLL in the K(K') val- 
leys have the four-component form = (<fio. m , 0, 0, 0) 
and ^o^ TO = (0, 0, 0), respectively, where 0o,m is 
the mth Landau orbital in the lowest (nonrelativistic) 
Landau level (see below). [By contrast, in the case 
of a real magnetic field, the second of these would be 
^am = (0, 0, 0, <f>o,m)-] Thus, in a pseudo-magnetic field, 
the PLL orbitals are counter-propagating, whereas in 
a real magnetic field they are co-propagating. Conse- 
quently, for a real magnetic field, the valley index can 
be regarded as a mere flavor degree of freedom, and a 
uniform valley-mixing perturbation mixes the Landau or- 
bitals in the two valleys. However, for strained graphene, 
the valley index should not be interpreted as a flavor in- 
dex, because the wavefunctions in the two valleys are 
spatially distinct. As we shall see, this implies that inho- 
mogeneous mass terms (i.e., mass terms with defects) are 
necessary to gap the LPLL. The precise consequences de- 
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pend on which gauge in pseudo- vector potentials is simu- 
lated by the strain pattern. (This is possible because, 
while two gauge-equivalent electromagnetic vector po- 
tentials arc physically identical, different pseudo-vector 
potentials are physically distinct as they correspond to 
different patterns of strain.) 

Uniform perturbations. We now address point (1) 
from the previous paragraph, ignoring spatial structure 
and considering the effects of uniform (i.e., defect-free) 
3CDW perturbations precisely at zero momentum (i.e., 
at the Dirac point). Many of the properties of the LPLL 
follow from the triviality of its sublattice structure: in 
particular, one can find the form of any perturbation in 
the LPLL by projecting it onto the A sublattice. Thus, 
perturbations within the LPLL are completely described 
by 2 x 2 matrices in r (i.e., valley) space (Table I). As 
a consequence, several perturbations that open up gaps 
in unstrained graphene, such as an intra-unit-cell charge- 
density wave ((X3T0) and the Kekule distortion {<J\T\), are 
trivial when projected to the LPLL. However, the SU{2) 
gauge potentials do open up gaps within the LPLL, 
as they project onto either t\ or r 2 (Table I), and can 
thus mix LPLL orbitals from the K and K' valleys. 

The two further perturbations that split the LPLL 
are valley polarization, m p = t^o-q, and the Haldanc 
mass [17], m# = T303. Within the LPLL, these terms 
are equivalent; both correspond to r 3 , which shifts the en- 
ergy of one valley with respect to the other. Both terms 
anticommute with the gauge fields; therefore, within the 
low-energy theory, it seems that these masses can be con- 
tinuously deformed into one another. (Thus, the fate 
of the topologically-protected edge mode associated with 
rriH cannot be addressed within the low-energy theory. 
Numerical calculations of the band structure in the pres- 
ence of both m,H and suggest that as Cx,j is increased, 
the edge state drifts away from the Dirac points; it is 
therefore plausible that £x,i and mu compete away from 
the Dirac points. We shall revisit this question in future 
work.) 

The perturbations discussed above are likely to be 
dynamically generated in experiments with neutral 
graphene (i.e., a half-filled LPLL), as opening up a gap 
would lower the ground-state energy. The 3CDWs could 
arise either because of the electron-phonon coupling or 
the electron-electron coupling; moreover, a spontaneous 
valley polarization might arise due to electron-electron 
interactions [8] . Alternatively, one can externally impose 
these gaps by growing graphene on an appropriately pat- 
terned substrate. 

Momentum dependence. We now turn to the second 
distinctive feature of the PLLs [point (2) above] and dis- 
cuss the role played by the spatial structure of PLL wave- 
functions. In order to address this, we shall consider the 
nature of the 3CDW perturbations on orbitals away from 
the Dirac point. 

We first consider the strain pattern realizing the 



Landau gauge; here, if the 3CDW amplitude is uni- 
form, the left- valley Landau orbital (po,m{x,y) oc 
gi^pz e -{y~ %J Tri ) H ^ being the magnetic length) 
can hybridize only with the right-valley Landau orbital 
cj)Q_ m (x,y) oc e % ^r x ^^^t' 2 ) /' 2 ; d ue to the conser- 
vation of k x . Although the overlap between these two 
states is nonzero, it decreases as |m| is increased, and 
becomes exponentially small for ml 3> L. This decrease 
of overlap leads to the convexity of the LPLL gap shown 
in Fig. 2, and implies that the 3CDW gap is an inherently 
mesoscopic phenomenon. (However, note that all experi- 
mental realizations of pseudo-magnetic fields in strained 
graphene involve mesoscopic systems.) 

Point defects. If the strain realizes a symmetric gauge 
pattern, (j>o im (r,9) cx e lmB r m e - r I 1 ) as i n the experi- 
ments of Rcf . [6] , the consequences are even more strik- 
ing. Here, the only allowed orbitals have m > 0, because 
otherwise <fio, m is not normalizable. As a result, no uni- 
form perturbation can gap out the m 7^ orbitals. Thus, 
in order to open up a gap, it is necessary to consider con- 
figurations in which the perturbations have defects, i.e., 
edge dislocations in the 3CDW case. Within the LPLL, 
the two independent 3CDW orders are represented by T\ 
and r 2 (Tabic I). The 3CDW can support edge disloca- 
tions [Fig. 1(c)], which take the form of vortex solutions 
within the LPLL theory: 



A(r) [n cos(n9) + r 2 sin(n#)] 



A(r)e inf> \ 
A(r)e~ ine J 
(5) 

where A(r) is a function that vanishes at r — and is 
constant for r 3> Im ■ Such vortex solutions carry angular 
momentum, so that they can lead to hybridization of 
different orbitals in the LPLL: 



J d9dr ($5 )m , 0) A(r) [n cos(n<9) + r 2 sin(n0)] 
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(6) 



for vorticity n, the Landau orbitals with m < n will 
be gapped. As a result, a larger vorticity can gap out 
more Landau orbitals and will thus lead to lower energy 
(although a larger vorticity might also cost greater elec- 
trostatic or elastic energy). In this case, one expects 
the defects to be dynamically generated along with the 
3CDW patterns. 

Effects in higher PLLs. We now turn to the effects of 
the aforementioned perturbations when projected onto 
higher PLLs. In the higher PLLs, the sublattice struc- 
ture is not trivial, so that the six masses are distinct 
perturbations. These perturbations fall into two classes, 
depending on whether their sublattice component is Co or 
cr 3 . Perturbations of the former class split all the PLLs. 
For the valley polarization term, this is obvious by inspec- 
tion: a polarization of the form t z moves all the PLLs in 
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the K up and all those at K' down, thus opening up a gap 
of size 2a between any pair of PLLs. Similarly, degener- 
ate perturbation theory shows that the perturbations 
mix \&|^ m and (assuming, for simplicity, that the 

strain pattern realizes the Landau gauge). By contrast, 
m H and do not mix \&Jf m and *^ m , and thus pre- 
serve the double degeneracy of higher PLLs. However, 
these perturbations do mix \&J£ m and m , thus shift- 
ing the energy of the nth PLL by a 2 /(2Aoy/n), where 
Ao is the cyclotron energy scale. 

These considerations might influence which set of per- 
turbations is in fact dynamically generated. In particu- 
lar, if one of the nonzero PLLs is half-filled, the favored 
perturbations are those of the a class. By contrast, 
at half- filling of the LPLL (i.e., for neutral graphene), 
the above argument suggests that the cr 3 class of per- 
turbations is preferred, as these lower the energy of the 
filled negative-energy PLLs. It is not clear, however, 
whether this energy saving is outweighed by changes in 
band structure away from the Dirac points. 

Experimental aspects. We close by touching upon var- 
ious experimental considerations. The 3CDW ordering 
described here can either be realized spontaneously in 
the presence of strain, or imposed externally. In the for- 
mer case, the pattern of ordering can be detected easily 
via scanning-tunneling microscopy (STM); this technique 
is commonly used to study stripe ordering, e.g., in high- 
temperature superconductors (see, e.g., Rcf. [18]). Alter- 
natively, one can study the 3CDW perturbations by ex- 
ternally imposing them. This is easiest to do in the case 
of engineered systems, such as molecular graphene [7], 
where a 3CDW pattern such as CyS can be imposed by 
hand on the triangular network of adsorbate molecules 
(or, alternatively, for hexagonal optical lattices [19]). 
However, it might also be possible to realize it using sys- 
tems of graphene grown on anisotropic substrates, which 
favor 3CDW formation. 
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